Texture control in a pseudospin Bose-Einstein condensate 
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We describe a wavefunction engineering approach to the formation of textures in a two-component nonrotated 
Bose-Einstein condensate. By controlling the phases of wavepackets that combine in a three-wave interference 
process, a ballistically-expanding regular lattice-texture is generated, in which the phases determine the com- 
ponent textures. A particular example is presented of a lattice-texture composed of half-quantum vortices and 
spin-2 textures. We demonstrate the lattice formation with numerical simulations of a viable experiment, iden- 
tifying the textures and relating their locations to a linear theory of wavepacket interference. 
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Topological spin textures arise in magnetic materials [1], 
in director fields of liquid crystals [2], in field theoretic mod- 
els of particles [3], and in models of the early universe [4]. 
Multicomponent Bose-Einstein condensates (BECs) may act 
as analogues of these and other condensed matter systems, 
enabling the study of phenomena that may be otherwise inac- 
cessible to experimental investigation. 

Bulk rotation of scalar (single-component) BECs provides 
one method for creating a regular lattice of quantized vortices, 
associated with mass currents. The additional spin freedom in 
rotating multicomponent BECs suggest the existence of re- 
lated spin lattice-textures, which are indeed observed [5]. 

In this Letter we describe a complementary method for the 
controlled production of a lattice-texture in a multicompo- 
nent BEC that forgoes bulk rotation of the condensate by in- 
stead exploiting interference to produce expanding lattices of 
singly-quantized vortices in any or all components. An ad- 
vantage of this method is that the textures comprising the mo- 
tif may be directly determined by the vortex-lattice alignment, 
which depends on the controllable wavepacket phases. Exper- 
imental production of isolated textures has been demonstrated 
within pseudospin-^, spin-1, and spin-2 BECs [6-8]. Tex- 
tures associated with small numbers of vortices in the lattice- 
texture could subsequently be isolated from the dynamically- 
expanding lattice by retrapping them optically. 

Nonrotated single-component BECs accommodate regular 
vortex lattices, created by a three- wave linear interference 
process [9]. In this scenario, three initially- separated BEC 
wavepackets expand and interfere, in an analogous process to 
a Young's three-pinhole interferometer [10]. When the initial 
wavepackets are arranged at the corners of an equilateral tri- 
angle, the resulting lattice has honeycomb symmetry and can 
be thought of as a dynamically-expanding Abrikosov lattice (a 
hexagonal lattice containing unit-charge vortices of one circu- 
lation), interleaved with a second similar lattice of vortices 
with opposite circulation. It is probable that such honeycomb 
vortex-antivortex (VA) lattices have been experimentally gen- 
erated in BECs, although this has neither been recognized nor 
verified directly — see [11, 12]. The initial wavepacket phases 
establish particular xy-translations of the resulting 2D VA lat- 
tices [9]. By engineering the initial wavepackets in a two- 
component BEC, we exploit this phase dependence to produce 



VA lattices within each component of a two-component BEC, 
aligning them to form the dynamically-expanding lattice tex- 
ture. In the example in Fig. 1, the position-dependent state 
is represented by local Bloch vectors, which project the state 
onto the surface of a unit-radius Bloch sphere (Fig. 1 inset). 

Although a rotating BEC is extended in three spatial di- 
mensions, its resident vortex lattice is primarily 2D in nature, 
which restricts the topology of any associated texture. The 
VA lattice is similarly 2D in nature, because the three initial 
wavepacket locations together define a plane. In both cases, 
axial absorption imaging captures the 2D structure by project- 
ing the density parallel to the vortex core nodal-lines. The 
simplest texture arises when a vortex in one component co- 
incides with a vortex-free region in the second component to 
create a half-quantum vortex (HQV) [13], which is a coun- 
terpart to the Alice string in particle physics [14]; its detec- 
tion in 3 He-A is sought as evidence of a spin triplet super- 
conducting state in Sr2Ru04 [15]. Another example, of rele- 
vance here, corresponds to a vortex in one component aligned 
with a vortex of opposite circulation in the second compo- 
nent. This texture may be thought of as a lower dimensional 
counterpart to the "spin vortex" that arises spontaneously in 
quenched 87 Rb spinor condensates [16]. Additional textures 
in 2D geometries include baby skyrmions [17], merons [18], 




FIG. 1. (color online) A planar lattice texture is created by the 
interference of a pseudospin-^ BEC initially separated into three 
pieces. The hexagonal lattice, visualized with Bloch vectors in this 
numerical simulation, has a motif containing three textures: one half- 
quantum vortex of each sign and a spin-2 texture. 



and planar spin textures [19], whose study has led to insights 
in reduced-dimensional superfluid or ferromagnetic systems, 
such as in the case of the quantum Hall effect [20] . We de- 
scribe a mechanism of lattice-texture formation that extends 
to three or more arbitrarily located spatially- separated spinor 
BEC pieces with arbitrary initial coherent phases. Although 
related to the Kibble-Zurek (KZ) mechanism [21], we focus 
on the vortex generation process in limits where linear inter- 
ference applies. A fuller description of the KZ mechanism 
might involve additional description of subsequent vortex dy- 
namics as part of an unspecified thermalization process. 

In the following, we present simulations of a proposed ex- 
periment to generate the lattice texture via the interference of 
a nonrotating pseudospin-^ BEC, initially divided into three 
pieces. Although it is possible to fully control the piece phases 
and resulting textures, we present a simplified experiment in 
which the phases are fixed. A planar hexagonal lattice-texture 
results, with a motif composed of three textures: one HQV 
of each sign and a spin-2 texture. We classify these textures 
by their topologies, and present a lattice model whose time- 
dependent growth is related to the initial BEC geometry. 

We numerically model a two-level 87 Rb BEC system with 
\F = l,m F = -1) = |1> and \F = 2,m F = +1) = |2> using a 
mean-field approach. These hyperfine states and their cou- 
pling behavior have been well- studied [22, 23] and are con- 
venient for the study of pseudospin-^ condensates. We per- 
formed 2D simulations, corresponding to pancake conden- 
sates. This geometry has the advantage that axial ballistic 
expansion is rapid, due to the initial tight axial confinement. 
Nonlinear effects, which might otherwise lead to distortion 
of any lattice and curvature or reconnection of vortex lines, 
are consequently short-lived. The order parameter field of 
a single-species BEC whose atoms occupy two internal hy- 
perfine levels is a position-dependent two-component pseu- 
dospinor ^(r) = (^(r),^^)), where r is a position vector 
[24] . The dynamical evolution of the BEC is governed by two 
coupled Gross-Pitaevskii Equations (GPEs): 
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where ¥* is the 2D order parameter of component |/) and 
V 2 is a 2D Laplacian. The self-interaction parameter Utj = 
4nti 2 aij/m depends on the intra- and inter-component s- wave 
scattering lengths a^ and the mass m of an atom of the con- 
densed species; there are three independent scattering lengths, 
since an = a<i\. A time-dependent factor r(0 results from the 
reduction from 3D to 2D, as described below. The normaliza- 
tion condition is £i J l^l 2 dv - N, where the total number of 
atoms N is preserved independently of the internal spin state. 
We use the scattering lengths reported by Mertes et al. [23]. 
The state-dependent scattering lengths are a\\ = 100.40 ao, 
ci22 = 95.00 ao, and 021 = an = 97.66 ao, where ao is the Bohr 
radius. The condensate contains N = 50, 000 atoms of 87 Rb, 
each of mass m = 1.4188 x 10" 25 kg. This relatively small 
population size was chosen to minimize nonlinear perturba- 
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FIG. 2. (color online) 87 Rb hyperfine levels |F = l,m F = -1) = |1> 
and \F - 2,m F - +1) = |2> are optically coupled with Gaussian and 
Laguerre-Gauss (LG) lasers, whose intensity profiles are shown. The 
LG^q phase is also shown winding from -n (black) to +n (white). 



tive effects and maximize the regularity of the lattice. 

Initially the BEC is tightly trapped in the axial direction, 
producing a pancake geometry, and further divided into three 
pieces by trapping within three transverse Gaussian potential 
wells, such as would be formed by three red-detuned lasers. 
The initial equilibrium condensate profile is established with 
all atoms in 1 1) by numerically evolving Eq. (1) through imag- 
inary time. In simulations, the numerical procedure estab- 
lishes a uniform phase (ip\ = (f2 = <f3) for the wavepack- 
ets, which fixes the translation of the resulting lattice texture. 
However, in typical experiments the initial wavepacket phases 
are random, resulting in the lattice being randomly translated. 
Therefore, in an experiment, both components must be im- 
aged simultaneously to correctly reconstruct the lattice texture 
structure. This might be performed by dual state imaging [25] 
or by nondestructive phase contrast imaging [26]. 

Having established the initial profile, we immediately re- 
move all trap potentials V; at t - 0. We account for ballistic 
expansion along the axial direction by assuming each piece to 
be a non-interacting gas with a Gaussian axial profile. Upon 
trap removal, the resulting rapid axial expansion and associ- 
ated rapid decrease of the nonlinear coupling term are mod- 
eled by the time-dependent factor T(t) = [mco/(2nfi)] l/2 (1 + 
t 2 co 2 )~ l/2 , where to is the angular frequency of the axial har- 
monic trap initially confining the condensate. We specify 
co = 500 rad s _1 , corresponding to a pancake geometry. 

We simulate the application of a two-photon n/2 pulse to 
excite half the atoms to |2> [Fig. 2]. This models an opti- 
cal process involving two lasers, coupling |1) and |2) via an 
intermediate level. The n/2 pulse may be applied either im- 
mediately before or after trap removal. The overall transla- 
tion of each VA lattice depends only on the relative phases of 
the initial wavepackets [9]. Although intra-component phases 
are uncontrolled in typical experiments, full control over the 
inter-component phase may be realized by spatially-localized 
Raman beam pairs focused on each wavepacket. Atoms ex- 
cited to 1 2) thereby acquire relative phases that produce a 
lattice which is predictably aligned with the lattice in |1), 
enabling the production of a continuum of related lattice- 
textures. In our simplified example, we instead employ a 
Laguerre-Gauss LG^Zq mode in one of the coupling beams 
to establish wavepacket phases cp\ + A<p - 27r/3, ^2 + A<p, and 
(f3 + A(f + 27r/3 for atoms in |2) [27]. The phase offset A<p 
has no effect on lattice translation [9]. The LG beam wave- 
front confers phase gradients on each wavepacket in |2). Any 
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FIG. 3. (color online) Numerical simulations of a two-component 
three-piece BEC generate vortex- antivortex lattices as the wavepack- 
ets expand and interfere. For the initial phases shown, the lattices 
align as shown in Fig. 4(a). 
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FIG. 5. (color online) Three texture types within the motif, (a) Vor- 
tex (v) and antivortex (av) cores are overlaid on the combined BEC 
components, (b) Bloch vectors are shown normalized to the com- 
bined probability density, (c) Order-parameter- space coverage for 
the three circular regions, where shading indicates either one or two 
overlapping projections. 



effects of this nonuniformity are minimized by virtue of the 
initial tight transverse confinement of the wavepackets. 

The BEC then evolves governed by Eq. (1), resulting in the 
axially-projected probability densities l^l 2 shown in Fig. 3. 
A honeycomb VA lattice is formed in each component, re- 
sulting from interference of the expanding wavepackets. The 
equal initial phases (cpi =(f2 = (f3) produce the particular lat- 
tice translations shown. The final frame of this figure shows 
the individual components l^l 2 at t = 43 ms, corresponding 
to the lattice texture in Fig. 1. 

We now consider the three textures that combine in the mo- 
tif of Fig. 1, corresponding to the different alignments of vor- 
tices shown schematically in Fig. 4(a). We examine these 
textures in more detail in Fig. 5, using simulation results at 
t = 69 ms. The locations of vortices and antivortices in each 
component are identified in Fig. 5(a). Because textures extend 
over a local neighborhood of the field, we specify three circu- 
lar bounded regions, labeled (i), (ii), and (iii), centered on the 
vortex cores. Each boundary is chosen to lie approximately 
halfway to the nearest neighboring vortex core. The same 
boundaries are shown in Fig. 5(b), which plots the field of 




FIG. 4. (a) Lattices containing vortices (white circles) and antivor- 
tices (black circles) form in components *Fi and ¥2 • Three combina- 
tions form different textures in *F = Q¥\ , ¥2). Gray circles indicate an 
aligned vortex and antivortex in the two components, (b) Wavepack- 
ets with centers given by ri = 0, r 2 , r 3 have relative phases 0, 02, 03, 
respectively. The resulting vortex- antivortex lattice basis vectors are 
r' 2 and r^, with a motif containing a vortex and antivortex separated 
by r' Q . Fringe spacing / is a convenient measure of unit cell size. 



Bloch vectors with lengths scaled by local density, and iden- 
tifies the three texture types. In crystallographic terms, the 
three textures constitute a motif, within a hexagonal lattice. 

In the two-component BEC system described here, the 
Bloch vectors within each bounded region of the physical 
space are projected onto a unit Bloch- sphere to reveal the 
texture topology. Maps of this order-parameter space are ob- 
tained by grouping adjacent unit Bloch vectors within the in- 
dicated circular regions, and coloring the patch of the sphere 
surface onto which they project. The maps are shown in 
Fig. 5(c) using a flattened equal-area sphere projection. 

Region (i) in Fig. 5 corresponds to a spin-2 texture with 
a zero net mass current, associated with the two constituent 
unit charge vortices of opposite sign that reside in different 
components. In traversing a small closed contour about the 
vortex core in Fig. 5(b) once, the vector projections wrap the 
Bloch sphere twice near the equator; the winding number is 
2. As contours of larger radius are traversed, the projections 
deviate from the equator toward the poles of the Bloch sphere. 
This texture is similar to spin textures described elsewhere in a 
system of 2D spins [19] and in three-component spin textures 
[16]. The two prominent holes visible in the coverage map 
[Fig. 5(c)(i)] are due to the influence of neighboring textures. 
These holes continue to shrink as the lattice expands. Cover- 
age of polar map-regions occurs when matter in one compo- 
nent coincides with a complete absence of matter in the other 
component, i.e., at vortex cores. Since the region (i) boundary 
is chosen to exclude neighboring vortex cores, the polar re- 
gions are not covered. The vortices associated with the spin-2 
texture are slightly misaligned, warping the texture. 

Regions (ii) and (iii) in Fig. 5 are centered on a vortex in 
one or other of the components. These correspond to HQVs. 
The sign of the HQV manifests as coverage of the correspond- 
ing hemisphere of the Bloch sphere [Fig. 5(c)(ii)— (iii)] . A 
threefold symmetry is apparent near the equator, due to the 
hexagonal symmetry of the lattice texture. As the lattice ex- 
pands, the projections on the boundary approach the equator. 
By identifying all points on the equator with a single value, the 
order-parameter space is restricted to one hemisphere. The 



HQVs are then topological defects, since the hemispherical 
order parameter spaces will be covered once by similar bound- 
aries to those shown. An HQV has wrapping number 1 for a 
hemispherical order parameter space, associated with its con- 
stituent unit charge vortex. It covers a 2n solid angle and is 
thus a 27r-defect whose sign depends on both the vortex sign 
and its resident component. Similarly, the boundary projec- 
tions of the spin-2 texture asymptotically approach lines of 
longitude on the Bloch sphere. If these linear boundaries are 
"healed," the wrapping number of the sphere is 2 and we may 
identify the texture as an 87r-defect. 

Having discussed the topology of the isolated textures, we 
now describe the lattice texture. The isolated textures expand 
along with the evolving lattice according to previously de- 
scribed expressions for vortex core locations in a model of 
interfering Gaussian wavepackets [9] . For vortices to be pro- 
duced by interference, the contributions from each expanding 
wavepacket must be approximately equal. This establishes a 
spatiotemporal condition, which limits the lattice extent and 
describes its growth by vortex formation within an expanding 
circular envelope [9] . By evaluating the locations of two ad- 
jacent vortices of the same sign within this envelope, we find 
that the lattice basis vectors [Fig. 4(b)] have lengths 



7i/(r 3 asm6), 



n/(r 2 a sin 6), 



(2) 



where r ; = |r 7 -| are the source spacings, and 6 is the interior 
angle at the origin of the triangle of wavepackets [Fig. 4(b)]. 
These expressions contain a time-dependent lattice scaling 
factor a = mht/[2(fit) 2 + 2m 2 (h/Ap) 4 ], which assumes that 
the wavepackets share a single initial momentum uncertainty 
A/?, and are therefore all of the same initial size. The distance 
/ in Fig. 4(b) corresponds to the bright fringe spacing, mea- 
sured experimentally by Henderson etal [12]. Forr = r 2 = r^ 
and = 7r/3, / = Ir^l V3/2 = n/(ra). The spacing of adjacent 
vortices and antivortices in the VA lattice is 

|ro| = 7i{r 2 2 + 2r 2 r 3 cos# + r 3 2 ) 1/2 / (3r 2 r 3 ar sin 0). (3) 

For r = r2 = r 3 and 6 = n/3, this reduces to Ir^l = 2n/(3ra). 

We have shown with numerical simulations that a lattice- 
texture forms in a nonrotated two-component BEC, initially 
separated into three pieces, and subsequently allowed to ex- 
pand and interfere. We created an expanding hexagonal 
lattice-texture by employing spatiotemporal wavefunction en- 
gineering to determine the specific relative phases of the ini- 
tial pieces and arrange them at the corners of an equilateral 
triangle. By employing a Laguerre-Gauss beam in the phase- 
engineering task, the lattice was created with a motif com- 
posed of two half-quantum vortices of opposite signs and one 
spin-2 texture, both being examples of direct relevance in re- 
lated condensed matter systems. More generally, the method 
presented provides a means for the deterministic production 
of a continuum of related lattice-textures. 

G.R. thanks T.P. Simula, K. Helmerson, L.D. Turner, R.P. 
Anderson, and E.J. Mueller for helpful discussions, and ac- 
knowledges Australian government support. 



* gary.ruben@monash.edu 
[1] S. Muhlbauer, B. Binz, F. Jonietz, C. Pfleiderer, A. Rosen, 

A. Neubauer, R. Georgii, and P. Boni, Science, 323, 915 

(2009). 
[2] C. Noel, Acta Polymerica, 48, 335 (1997). 
[3] T. H. R. Sky rme, Proc. R. Soc. A, 260, 127 (1961). 
[4] A. Vilenkin, |phys. Rev. D, 23, 852 (1981) . 
[5] T. Mizushima, N. Kobayashi, and K. Machida, P hys. Rev. A, 

70, 043613 (2004); J. W. Reijnders, F. J. M. van Lankvelt, 

K. Schoutens, and N. Read, ibid., 69, 023612 (2004); E. J. 

Mueller, ibid., 69, 033606 (2004). 
[6] M. R. Matthews, B. P. Ander son, P. C. Halian, D. S. Hall C. E . 



Wieman, and E. A. Cornell, Phys. Rev. Lett., 83, 2498 (1999) 



[7] A. E. Leanh ardt, Y. Shin, D. Kielpinski, D. E. P ritchard, and 



W. Ketterle, Phys. Rev. Lett., 90, 140403 (2003) 



[8] L. S. Leslie, A. Hansen, K. C. Wright, B. M. Deutsch, and N. P. 

Bigelow, Phys. Rev. Lett., 103, 250401 (2009). 
[9] G. Ruben, D. M. Paganin, and M. J. Morgan, Phys. Rev. A, 78, 

013631 (2008). 
[10] G. Ruben and D. M. Paganin, frhys. Rev. E, 75, 066613 (2007J ; 



76, 029901(E) (2007) 



[11] D. R. Scherer, C. N. Weiler, T. W. Ne ely, and B. P. Anderson, 

Phys. Rev. Lett., 98, 110402 (2007)| ; R. Carretero-Gonzalez, 

B. P. Anderson, P. G. Kevrekidis, D. J. Frantzeskakis, and C. N. 

Weiler, Phys. Rev. A, 77, 033625 (2008). 
[12] K. Henderson, C. Ryu, C. MacCormick, and M. G. Boshier, 

New J. Phys., 11, 043030 (2009). 
[13] M. C. Cross and W F. Brinkman, J. Low Temp. Phys., 27, 683 

(1977). 
[14] U. Leonhardt and G. E. Volovik, JETP Lett., 72, 46 (2000); 

T. Isoshima, K. Machida, and T. Ohmi, J. Phys. Soc. Jpn., 70, 

1604 (2001). 
[15] V. Va karyuk and A. J. Leggett, P hys. Rev. Lett., 103, 057003 

(2009). 
[16] L. E. Sadler, J. M. Higbie, S. R. Leslie, M. Vengalattore, 

and D. M. Stamper-Kurn, Nature, 443, 312 (2006); H. Saito, 

Y. Kawaguchi, and M. Ueda, Phys. Rev. Lett., 96, 065302 

(2006). 
[17] B. M. A. G. Piette, B. J. Schoers, and W J. Zakrzewski, Z. 

Phys. C, 65, 165 (1995). 
[18] C. G. Callan, R. Dashen, and D. J. Gross, Phys. Lett. B, 66, 

375 (1977). 
[19] N. D. Mermin, Rev. Mod. Phys., 51, 591 (1979). 
[20] T. Ando, Y Matsumoto, and Y Uemura, J. Phys. Soc. Jpn., 39, 

279 (1975). 
[21] T. W B. Kibble, J. Phys. A, 9, 1387 (1976). 
[22] M. R. Matthews, D. S. Hall, D. S. Jin, J. R. Ensher, C. E. Wie- 
man, E. A. Cornell, F. Dalfovo, C. Minniti, and S. Stringari, 

Phys. Rev. Lett., 81, 243 (1998). 
[23] K. M. Mertes, J. W Merrill, R. Carretero-Gonzalez, D. J. 

Frantzeskakis, P. G. Kevrekidis, and D. S. Hall, Phys. Rev. 

Lett., 99, 1 90402 (2007). 



[24] T.-L. Ho, E>hys. Rev. Lett., 81, 742 (1998) 



T. Ohmi and 

K. Machida, J. Phys. Soc. Jpn., 67, 1822 (1998). 
[25] R. P. Anderson, C. Tickno r, A. I. Sidorov, and B. V. Hall, P hys. 

Rev. A, 80, 023603 (2009). 
[26] L. D. Turner, K. F. E. M. Domen, and R. E. Scholten, Phys. 

Rev. A, 72, 031403(R) (2005). 
[27] M. F. Andersen, C. Ryu, P. Clade, V. Natarajan, A. Vaziri, 

K. Hel merson, and W D. Phillips, P hys. Rev. Lett., 97, 170406 

(2006) = 



